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We present an approximative simulation method for quantum many-body systems based on coarse 
graining the space of the momentum transferred between interacting particles, which leads to effec- 
tive Hamiltonians of reduced size with the flavor-twisted boundary condition. A rapid, accurate, 
and fast convergent computation of the ground-state energy is demonstrated on the spin-^ quantum 
antiferromagnet of any dimension by employing only two sites. The method is expected to be useful 
for future simulations and quick estimates on other strongly correlated systems. 
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Understanding a quantum many-body system is fundamentally challenging because of the exponential growth of the 
number of states with the system size. To estimate the physical properties of an intractably large system, a common 
practice is to extrapolate the results for several significantly reduced sizes^ together with the periodic boundary 
condition (PBC)i^ A different yet complementary approach is to continuously twist the boundary conditions for one 
solvable size. The typical flavor-independent version of the latter has been seen in solid state physics,^ while the 
flavor version has been used in quantum chromodynamicsi^'i Here we apply the flavor-twisted boundary conditions 
(FTBC)^^ to the spin-i quantum antiferromagnet, one of the basic models in solid state physics. We show that the 
ground-state energy can be accurately calculated using only two sites. 

We begin with an explicit derivation of FTBC, which is necessary for systematical studies of this and related 
methods. Let us consider the connection between a large system of size L and a small one of size I, both with 
PBC. (To distinguish them, the upper-case and lower-case letters are used thoroughly for the large and small lattices, 
respectively.) The large system is described by the following general Hamiltonian in the second quantized language 
and in the real space representation, 

Hl = J2 ^ClC 3a + h.c.) + £ U& n/ ClC J/a ,C I1(r ,C Ja , (1) 

I,J,cr I.J.a 

j'jV 

where C Icr denote the quantum operator that annihilates a particle with flavor a at site I. L = Yl a =i ^ a ^ s ^ ne 
total number of the lattice sites in a d-dimensional space with L a being the site number in the a-dimension. The 
thermodynamic limit is reached when all L a — ► oo. 

There also exists a reciprocal space where the counterpart of the site I is the momentum point K. Imposing PBC 
to the real space discretizes the K space as follows^ 

K a = j^-M, M = 0, 1, 2, . . . , L a - 1 (2) 

where all the lattice constants have been scaled to unit. Eq. ^ translates the concepts of large and small in the real 
space to those of fine and coarse in the K grid, respectively. In the reciprocal space, 

Hl = E^^Ko-^Kct + ^- E ^k,k',q ( ^k+q, ( tC'k'-q.<t' ( ^k',o''C'k,o-} ( 3 ) 

Kcr Q Kcr.K'cr' 

where C K(T annihilates a particle with momentum K and spin a and is given by 

Ck. = ^ E e " K ' RlC i- ( 4 ) 
where Ri is the coordinates of the I-th site. The bare energy dispersion and the interaction function are 



= E *&e*-C*-*>. (5) 
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FIG. 1: Illustration of the momentum transfer due to the particle-particle interaction: (a) original and (b) coarsened. 



Eq. §3§ appears to describe two particles with K and K', respectively, interact with internal momentum transfer 
Q, as diagrammed in Fig. [TJa). Since the Q points are to be integrated [c.f. X^q{' ' ' } m Eq. ©], here comes a well- 
known numerical trick: the summation (integration) over a fine grid may be well approximated by that over a rather 
coarse grid. This numerical recipe receives particular support in the quantum many-body systems of interest, where 
a small number of Q are far more important than the others, the so-called Q-mode resonance, e.g., Q = (ir, ir, ■ ■ ■ , ir) 
in antiferromagnetic correlation. The Q-mode resonance implies dramatic response to small stimuli such as changing 
temperature, applied fields, pressure, doping, etc, which is important to determining the functionalities of the system. 
Therefore, it may be a good approximation to coarsen the Q grid as long as the important Q points are included. 

To find a way to coarsen the Q grid such that the resulting Hamiltonian can be readily transformed back to the 
real space (where the parameters are often short-ranged, meriting for numerical computation), we consider the small 
system of size / with PBC whose momenta k are given by: 

k a = — m, m = 0, 1, 2, . . . , l a — 1. (7) 

i a 

where I = n^ =1 Z a <C L is the total site number of the smaller lattice. Let the small system be commensurate with 
the large one, that is, all the k points can be found in the K grid of the large system. Then, any momentum in the 
fine K grid can be rewritten as, 

K = K + k, K' = K' + k', Q = Q + q, (8) 

where k' and q arc in the k grid defined in Eq. ([7]); the L/l 'twists' K (likewise K' or Q) consist of a subset of the 
original K points that fill the interstitial region of the k points nearby the origin. For example, for d = 1, 

K a e(-f,f}. (9) 

As another example, I = 2 for the square lattice means that k = (0, 0) and (ir, tt), and that K cover half of the original 
K space satisfying cos^a,) + cos(K y ) > 0. 

One can explicitly rewrite Hl in terms of the real-space Z-site lattice by introducing a new set of particle annihilation 
operators, 

d k .;K = C K +k . = C K CT , (10) 

'^k ^y.'' kr <iu,.K- (n) 

where is the coordinates of the i-th site of this /-site lattice, satisfying 



E 
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= lS i; y (12) 
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By using Eqs. ([7l fl2|) , the original Hamiltonian for the L-site system, Eq. ([I) or Eq. ([3]), can be transformed to a set 
of /-site subsystems described by the following new Hamiltonian in the real space representation 

Ker >J Q Kcr,K'<r' i.i'jj' 

f e iK-(r 3 - ri ) «'.(r 4 ,-.y) jQi'i'-^U^'d] - ,d\, -, -d., lici d- it)- (13) 
t U J 1 i<r;K+Q jV':K'-Q i'cr';K' j<t;K J > 

This establishes the exact transformation between the large and the small, subject to the linear size of the small 
system being not shorter than the range of any Hamiltonian parameters of Eq. (fTj). 
Approximation 1. — Now we coarse the Q grid with 

Q = 0. (14) 

This means that only the q points, whose total number is I, are retained as the internal momentum transfer. Then, 
K + q = K + k + q = K + k'. That is, any K point is scattered to the sum of the same twist and another point in 
{q}. Thus, a particle with a given K can be scattered by the particle-particle interaction into only I points in the K 
momentum space, instead of L points, as illustrated in Fig. [ljb). Eq. (|13|) is rewritten as 

Hl ~1 E Ht(Kv,KW), (15) 

Ka.K'a' 

where the subsystems are 

fli(K4V) = Et elIC ' (r3 " r,)< ij rf L ; K d j.;K + /i - C -] + 



Z-^i v y icr':K' jo-';K' ' 



jo 
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+ y e^-^e^^'-^U^d] -4 , - d., !it! d. ft- (16) 

/ j 'J J 1 ict;K jV;K' i'«r;K' jer;K v ' 

UJ'.i' 

Combined with Eq. ()15p . Eq. Q16p for each pair of K and K' is not independent of another pair. For example, the 
{K, K'} subsystem and the {K, K"} subsystem share the same momentum points with K. Therefore, all the (L/l) 2 
subsystems of size I are connected and one still has to deal with a case of size L. To simplify the calculations to size 
I, a further approximation is needed. 

Approximation 2 (FTBC). — A simple approximation is to treat iJ;(Kcr, K'er') independently. Thus, to estimate the 
expectation value of an observable Ol in the large system is to calculate 

traced 6 L ) ~ ^ traceO^.^^' A;K<t,k' ct ')> (17) 

Ko-,K'cr' 

where ph denotes the density matrix for the large system and P;.^ ^, a , denotes the density matrix for the isolated 
small system with the twists K and K', as determined by Eq. (HHJ). kV ^ s * ne transformed observable following 

Eqs. (fl"Tj|) - (fl"2]) and Eq. (fM)) . For example, the ground state energy is approximated by 

^(Ka.KV)!^)- £ (^^jH^K'a')^^,), (18) 

Kct.K'ct' Ker,K'er' 

where |$^) denotes the ground-state wave function for the large system and |</>;.k<t kv) denotes a wave function for 
the isolated small system with the twists K and K'. 

This approximation is actually equivalent to FTBC (Ref. 0) as explained below. Comparing Eq. p6[) with Eq. (JTJ) , 
one finds that they are very similar in form — only differ in the parameters by a phase factor e' K '( ri-r j) associated 
with the momentum twist K. Actually, solving Eq. (|16p is equivalent to solving the following Hamiltonian for the 
i-site lattice 

H((Ktr,KV) = V(tf,ttf #d. ^ + h.c.)+Y [t?-d] ,~d. ,^,+h.c.) 

>V ' ' V IJ i CT; K jer;K ' v 'J ict';K' J<t ';K' ' 

ij ij 

+ V Uwwd *4 ,*A ,*,<*. ,> (19) 
Z_/ 11 JJ icr:K i'er':K' jV;K' j<t;K ' 

ij j'.i' 
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together with the boundary conditions such that translating its many-body wavefunction |\I>;) along the a-th dimension 
l a steps yield e^ NlKa+N2K ^ la \H>i) where Ni and N? are the numbers of the two flavors of particles associated with 
the two twists K and K', respectively. 

With FTBC, the problem of solving the original Hamiltonian for a L-site system is reduced to solving (L/l) 2 
Hamiltonians for /-site subsystems, each corresponding to a given pair of K and K'. The computational benefit 
is substantial, since the number of states grows exponentially with the system size. For the Hubbard model as an 
example, the calculation load is reduced from 0(4 L A L ) to 0[(L/l) 2 x 4 Z 4']. In addition, the (L/l) 2 prefactor is fully 
parallelizable and it can be readily handled by using the same integration trick of replacing a fine-grid sum with 
a coarse-grid sum. Since FTBC is reached at the level of Hamiltonian and interpreted as boundary conditions, it 
has least limitation and fully compatible with other many-body approaches, e.g., Lanczos exact diagonalization and 
quantum Monte Carlo methods^ 

Approximation 3. — The typical twisted boundary conditions^ (TBC) can be reached from Eqs. (fT5")) and (TIT))) by 
further taking the K = K' approximation, 

i ]T #,(K<7,KV)« H^Ko-^o-'). (20) 

Ker.K'er' Ka.a' 

With TBC, each subsystem has only one momentum twist K. Commonly, the one-site or two-site (ry = and 
Tji = Tj) interactions are the leading interaction terms. Then, the twisted phases of the interaction terms in Eq. (p~6|) 
cancel, which renders the interaction terms for I to be of the same form as those for L. Thus, the approximation by 
TBC allows the continuous sampling of the momentum space for one-particle excitations, but it prohibits the same 
for two-particle excitations. To compare, FTBC allows both in principle. It could be expected that FTBC is more 
accurate, as illustrated below. 

To complete, PBC is an additional approximation of TBC (K = K' = 0). FTBC with one twist zero fixed (K' = 0, 
referred to as FTBC0)££ was studied before. 

To illustrate all the above points, we use FTBC to estimate the ground state energy of the spin-1/2 Heisenbcrg 
quantum antifcrromagnct of any dimension. Here the flavor is the spin of electrons, consisting of | and J..— The results 
are compared with those obtained from using the other boundary conditions and linear spin- wave theory (LSW).— 

The spin-1/2 antiferromagnetic Heisenberg model is given by 

h l = jJ2 Pf s J + 1( s i s j + s j s j )]> ( 21 ) 
<I-J> 

where the spin operator Sj'~' z = J^^u ^i/i^/tr/ - ^iv w ^ n u ^lv~' z being the Pauli matrix elements. j) runs over 
nearest neighbors. It is not only the basic account of antiferromagnetismi^ but also the ground zero of understanding 
high-temperature superconductivity in copper oxidesi 11 ' 12 Two interesting states have been intensively considered in 
the literature: the long-range ordered Neel state and the resonating-valence-bond (RVB) state i 12 ' 13 The concept of 
RVB is based on the fact that the minimum bond energy (— 0.75J) is realized in a two-site system, much lower than 
the bond energy (—0.25 J) of the classic Neel state; the valence-bond between two nearest-neighbor spins is arguably 
key to understanding low-dimensional quantum antifcrromagnets J^. However, the 2-site dimer breaks bonding to spins 
on the other sites in an extended system. For example, the total energy from the dimers is — 0.375J per bond for 
d = 1 and —0.1875 J per bond for d = 2, much higher than the exact resul t 14 ' 15 — 0.443J for d = 1 and the numerical 
result^ —0.334 J for d = 2, respectively. It is argued that the resonating (i.e., superposition of the degenerate states 
of different dimer configurations) could lower the energy substantially^ But an accurate estimation of the ground 
state energy was achieved only when long-distance spin dimers were also included* 1 ^ 
With FTBC, after coarsening the Q mesh, one obtains 

Hi (Kf,Ki) = J^>f S J + he^-^y^-^s+sr + h.c.]}, (22) 
<iJ> 

where s? = d^^d.^, s7 = df^d,^ and sf = |(4 T; KT d iT;KT ~ rf U ; Ki d U;Kj.)- The spin-exchange terms strongly 
depend on the double twists, K| and K^, for the spin-up and spin-down electrons, respectively. In comparison, the 
result with TBC is 



(23) 
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TABLE I: The ground-state energy per bond (in unit of J) of the spin-1/2 Heisenberg quantum antiferromagnet in the 
thermodynamic limit for the linear chain (d = 1), the square lattice (d — 2), and the body-centered-cubic lattice (d — 3). The 
results from using FTBC, FTBCO, TBC or PBC are calculated with I = 2. 
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independent of twists, the same as with PBC in this case. 

The accuracy and convergency of FTBC are tested with the bond energy 
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"-bond \LJ J—t 
K T K i 

where ribond is the number of the nearest-neighbor AF bonds in the small system. Let us first take the most dramatic 
approximation, viz: I = 2 (note ribond = 1), the fundamental of the RVB state. With PBC, the q mesh contains only 
two points: (0, 0, ... , 0) and (n, w, . . . , n) corresponding to a spin singlet and a triplet, respectively; the energy of the 
q = state is —0.75 J. By using FTBC to continuously and smoothly twist the energies of the subsystems, the bond 
energy becomes 

E(L;l = 2) = J2 [-i-ie*<*T-ftiH«-*0] 



Kt -K i 



j j A i 
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{ II jhf\ dR tJ d^fcos(i^/2)cos(iq72) 



L (2tt) 2 _ 

for L — > oo. (25) 
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The results are listed in Table Q] together with those obtained from using FTBCO, TBC, and PBC as well as the 
LSW results. The bond energy in the thermodynamic limit and for any dimension is accurately reproduced using the 
simple Eq. (J2SJ for I = 2. 

Next, the convergency of the bond energy E(L — * oo;l) with respect to I is presented in Fig. [Ha). The energy 
of each subsystem is calculated with the Lanczos exact diagonalization method and a 5-point Gaussian integration 
in the twist momentum space, which is enough for the 6-digit accuracy. It is obvious that the estimate from FTBC 
converges to the exact solution — J In 2 + J/4 (Refs. [TH andll5[) much faster than the extrapolation to L — > oo from 
finite-size scaling of the PBC or TBC results. Finally, E(L; 2) and E(L\ 4) for a number of L are plotted in Fig. HJb) 
and compared with the exact solutions for L in order to show the fast convergence of FTBC with respect to the size 
of the simulated system. Overall, the errors for I = 4 are smaller than for / = 2 (0.003 versus 0.009 for L = 24). For 
I = 2, the small deviations from the exact solutions for L > 4 follow a power law -0.009510(5) + 0.152(2)//> 976 ( 6 ) 
with x 2 = 1-7326 x 10~ n . This means that the error grows rather slowly as L increases away from I. These results 
indicate that a large-scale feature of a quantum antiferromagnet could be captured at the length scale of one lattice 
constant with FTBC. 

Finally, it is worth mentioning that in the above derivation of FTBC, we have revealed a less approximated approach, 
i.e., Eq. (|14[) alone. It is interesting to explore this approach and compare it with FTBC. Also, the explicit formulation 
of FTBC could facilitate to devise other approximations with lighter computational load. These studies are beyond 
the scope of the present work and will be published elsewhere. 

Summarizing, based on coarse graining the space of the momentum transferred between interacting particles, 
we have derived the flavor-twisted boundary condition for simulating quantum many-body systems with effective 
Hamiltonians of reduced size. A rapid, accurate, and fast convergent computation of the ground-state energy is 
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FIG. 2: (a) Bond energy of the one-dimensional spin-1/2 Heisenberg antiferromagnet for L — > oo estimated from using FTBC 
(diamonds) and TBC or PBC (circles) for a number of I. The horizontal solid line denotes the exact solution ( — Jln2 + J/4). 
The dotted lines are guide to eyes, (b) Bond energy of the one-dimensional spin-1/2 Heisenberg antiferromagnet for a number 
of L estimated from using FTBC for I — 2 and 4 (triangles), respectively, compared with the exact solutions (circles). 



demonstrated on the spin-^ quantum antiferromagnet of any dimension by employing only two sites. The method is 
expected to be useful for future simulations and quick estimates on other strongly correlated systems. 
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